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Abstract 

We calculate the electric dipole moment of electron using the experimental result of 
muon electric dipole moment and upper limit of the BR{^ — > erf) in the framework of 
the general two Higgs doublet model. Our prediction is 10~ 32 e — cm, which lies in the 
experimental current limits. Further, we obtain constraints for the Yukawa couplings 
|jv re an d tu- Fhrally we present an expression which connects the BR(t — > fij) and 
the electric dipole moment of r-lepton and study the relation between these physical 
quantities. 
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1 Introduction 



Lepton Flavor Violating (LFV) interactions are one of the important source of physics beyond 
the Standard model (SM) and have reached great interest with the improvement of experimental 
measurements. The processes \x — > c~y and r — > /ry are the examples of LFV interactions and the 
current limits for their branching ratios {BR) are 1.2 x 10~ n |T[ and 1.1 x 10~ 6 respectively. 
In such decays, the assumption of the non-existence of Cabibbo-Kobayashi-Maskawa (CKM) 
type matrix in the leptonic sector forbids the charged Flavor Changing (FC) interactions and 
therefore the physics beyond the SM plays the main role, where the general two Higgs doublet 
model (2HDM), so called model III, is one of the candidate. In this model, LFV interactions 
can exist at loop level with the help of the internal neutral higgs bosons ho and A . The Yukawa 
couplings appearing in these loops are free parameters and their strength can be detemined 
by the experimental data. In the literature, there are several studies on LFV interactions 
in different models. Such interactions are studied in a model independent way in ||, in the 
framework of model III 2HDM ||, in supersymmetric models [|, §, 0, §, 0, [II], [TT|] . 

CP violating effects also provide comprehensive information in the determination of free 
parameters of the various theoretical models. In extensions of the SM, new source of CP- 
violating phases occur. The exchanges of neutral Higgs bosons in the 2HDM \T2] or the charged 



sector in multi Higgs doublet models (more than two Higgs doublets) [13fl induce CP violating 
effects. In the model III 2HDM, one of the source for the CP violation is the complex Yukawa 
coupling. Non-zero Electric Dipole Moments (EDM) of the elementary particles are the sign 
of such violation. These are interesting from the experimental point of view since there are 
improvements in the experimental limits of charged lepton EDM. EDM of electron, muon 
and tau have been measured and the present limits are d e = 1.8 ± 1.2 ± 1.0 x 10~ 27 [114 



dp = (3.7 ± 3.4) x 1(T 19 |15l and d T = 3.1 X 10 



-16 



Dipole moments of leptons have been studied in the literature extensively. In [[l?]] , it is 
emphasized that the dominant contribution to the EDM of lighter leptons comes from two loop 
diagrams that involve one power of the Higgs Yukawa couplings. In this work, the CP-violation 
is assumed to be mediated by neutral Higgs scalars and the EDM of electron is predicted 
at the order of the magnitude of 10~ 26 e — cm. Further EDM of leptons have been analyzed 
in supersymmetric models || [TT|, |T5| . In the recent work [2T1] EDM of leptons are studied 



by scaling them with corresponding lepton masses and the EDM of electron is predicted as 
10~ 27 e — cm. 

In our work, we study EDM of leptons e, fi, r and LFV processes \i — > ey, r — > fiy in the 
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general 2HDM (model III). The source of EDM of a particle is the CP violated interaction 
and it can come from the complex Yukawa couplings. In the model III, it is possible to get a 
considerable EDM at one-loop level. Further, LFV processes can exist also at one-loop level 
with internal mediating neutral particles h and A since there is no CKM type matrix and 
therefore no charged FC interaction in the leptonic sector according to our assumption. 

The paper is organized as follows: In Section 2, we present EDMs of e, /i, r leptons and the 
expressions for the decay width of processes \i — » and r — > /ry in the framework of the model 
III. Section 3 is devoted to discussion and our conclusions. 

2 Electric dipole moments of charged leptons and LFV 
interactions in the general two Higgs Doublet model. 

The 2HDM type III permits flavor changing neutral currents (FCNC) at tree level and makes 
CP violating interactions possible with the choice of complex Yukawa couplings. The Yukawa 
interaction for the leptonic sector in the model III is 

C Y = vSliL<l>iE j R + g j liL(f>2E jR + h.c. , (1) 

where i, j are family indices of leptons, L and R denote chiral projections L(R) = 1/2(1=F7 5 ), 0, 
for % = 1,2, are the two scalar doublets, l iL and E jR are lepton doublets and singlets respectively. 
Here 0i and 02 are chosen as 



^ = 7=2 



\ , ( V2 X 
v + H° ) + 1 i X ° 



and the vacuum expectation values are 



<^>=^[° v ) ;<0 2 >=o. (3) 

With this choice, the SM particles can be collected in the first doublet and the new particles 
in the second one. The part which produce FCNC at tree level is 

£y,fc = C7/./02 /•-',/,> + h.c. . (4) 

Here the Yukawa matrices have in general complex entries. Note that in the following we 
replace £ D with where "N" denotes the word "neutral". The complex Yukawa couplings are 
the source of CP violation and EDM of particles are created by the CP violated interactions. 
The effective EDM interaction for a lepton is defined as 

£edm = id t / 7 5 <V I , (5) 



where F^ v is the electromagnetic field tensor and "d" is EDM of the lepton. Here, "d" is a real 
number by hermiticity. Concentrating only on neutral currents in the model III, the neutral 
Higgs bosons ho, A Q can induce CP violating interactions which can create EDM at loop level. 
Note that, we take H 1 and H 2 in eq. (|2|) as the mass eigenstates h and A respectively since 
no mixing between CP-even neutral Higgs bosons ho and the SM one, Ho, occurs at tree level. 
Due to possible small mixing at loop level we also study the mixing effects on the physical 
parameters in the Discussion part. 

Now, we give the necessary 1-loop diagrams due to neutral Higgs particles in Fig. |I[ Since, 
in the on-shell renormalization scheme, the self energy J2{p) can be written as 

J2(p) = (p- m l )J2(p)(P ~ mi) , (6) 

diagrams a, b in Fig. [I] vanish when /-lepton is on-shell. However the vertex diagram c in Fig. 
|T] gives non-zero contribution. The most general Lorentz-invariant form of the coupling of a 
charged lepton to a photon of 4-momentum q v can be written as 

r„ = Gi(q 2 ) 7 M + +G 2 (q 2 ) (T/jj, q v 

+ G 3 (q 2 )a^ 5 q» (7) 

where q v is photon 4- vector and q 2 dependent form factors G\(q 2 ) and G2{q 2 ) are proportional to 
the charge and anamolous magnetic moment of /-lepton respectively. Non-zero value of G 3 (q 2 ) 
is responsible for the CP violation and it is proportional to EDM of /-lepton. By extracting 
this part of the vertex, /-lepton EDM "<i/" (/ = e, /i, r) (see eq.(^)) can be calculated as a sum 
of contributions coming from neutral Higgs bosons ho and A , 

di = -^f^ % &)* ~ mVHo) - F 1 (y Ao )) , (8) 

for I = e,\i and 

iGf e r Q 



^ = -^^{^((^; T ) 2 -(0 2 )(^(^o)-^2(^ )) 

Q, ^m:,r) 2 - (O 2 ) K In (z h0 ) - r Ao In (z Ao ))} , (9) 



where the functions Fi(w), F 2 {w) and F 3 (w) are 

w (3 — 4 w + w 2 + 2 In w) 



FAw) 



[-1 + w 

2{-2 + w)wln\{<Jw - y/w - 4) 



F 2 (w) = wlnw + — 2 = -. (10) 

Jw {w — 4) 



Here yn = 7^27, ru = 7^77 and zh = 7^727, Q T and Q M are charges of r and /x leptons respectively. 



In eqs. (|§|) and (g) is defined as £^ = y^^^N,ij- I n ec f © we take info account 
only internal r-lepton contribution since, in our assumption, the Yukawa couplings C^r/, h3 = 
e, [A, are small compared to ^ ri i = e,fi,r due to the possible proportionality of the Yukawa 
couplings to the masses of leptons underconsideration in the vertices. In eq. (|9|) we take also 
the internal /i-lepton contribution besides internal r-lepton one and ignore the one coming from 
the internal e-lepton respecting our assumption (see Discussion part). Note that, we make our 
calculations in arbitrary q 2 and take q 2 = at the end. 
Using the parametrization 

S.* = ISU^ . (ii) 

the Yukawa factors in eqs. (|8|) and (||) can be written as 

m; lT r - (O 2 ) = - 2 « 2 ^ \^j 2 (12) 

where I = e, /i, r. Here 81 are CP violating parameters which are the sources of the lepton 
EDM. 

Now, let us consider the lepton number violating process fi — > ej which is a good candidate 
in the determination of the Yukawa couplings and new physics beyond the SM. Since we take 
into account only the neutral Higgs contributions in the lepton sector of the model III, the 
contribution comes from neutral Higgs bosons ho and Aq (see Fig. [I]). In the on-shell renor- 
malization scheme the self energy diagrams are cancelled and the vertex diagram (Fig. |l|-c) 
contributes. Taking only r lepton for the internal line, the decay width Y reads as 

r( Ai ^ e7 ) = c 1 (|A 1 | 2 + |A 2 | 2 ) , (13) 



where 



At = Q T - Ce^( F l(^o)-^l(^o)) , 

A2 = Qt q ~ £>N*er ^N*fir (-^1 {l/ho ) -^1 {VAq ) ) > (14) 



C\ = F 3 2™4 M and the function F 1 (w) is given in eq. fllDfl. Here the amplitudes A\ and A 2 



have right and left chirality respectively. In eq. (|13|) we ignore the contributions coming from 
internal /i and e leptons respecting our assumption on the Yukawa couplings (see Discussion). 

Another LFV process is r —>■ /ry and it is rich from the theoretical point of view. The decay 
width of this process can be calculated using the same procedure and reads as 

r( r ^ Ai7 )= C2 (| jBl | 2 + | J B 2 | 2 ) , (is) 
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where 

Bl = Q' A*^ m ^{^(gl(^o)+^l(^ ))+6^, rT (Fl(y feo )-^l(^o)} , 

B * = Q t az1 m ^T{t%AGi{y hQ ) + G 1 (y Ao )) + 6^* rT (iq(^ ) - i^ftO} , (16) 

and C2 = F 327^ mT • Here the amplitudes £>i and 5 2 have right and left chirality respectively. 
The function G\{w) is given by 

w (2 + 3w — 6w 2 + w 3 + 6wlnw) 

(-1 + Wj 4 

Note that, in the case of the possible loop induced mixing of the neutral Higgs bosons ho and 
the SM one H , the functions Fi(y ho ), F 2 (r ho ), ln(z ho ) and G\{%)h ) in eqs. (|| [|, [14|, |16|) should 
be multiplied by cos 2 a where a is the small mixing angle. 

3 Discussion 

The Yukawa couplings i, j — e,fi,r play the main role in the determination of the lepton 
EDM and the physical quantities for LFV interactions, in the model III. These couplings are 
complex in general and they are free parameters which can be fixed by present and forthcoming 
experiments. The experimental results on the EDM of leptons e, fi, r and LFV process fi — > ej 
are our starting point of the predictions on Yukawa couplings related with leptons. Note that, 
in our predictons, we assume that the Yukawa couplings i, j = e,/i, are small compared 
to £jv T j« = e, fi, t since the strength of these couplings are related with the masses of leptons 
denoted by the indices of them, similar to the Cheng-Sher scenerio PTJ . Further, we assume 



that is symmetric with respect to the indices i and j. 

Since non-zero EDM can be obtained in case of complex couplings, there exist a CP violating 
parameter Q\ coming from the parametrization eq. (|TTD. At this stage, we find a restriction for 
%n,th usrn g ec L- © an d the experimental result of fi EDM p^j , 



0.3 x 1(T 19 e - cm < d M < 7.1 x 10~ 19 e - cm . (17) 

Fig. || shows mA dependence of £jy )T „ for sinO^ = 0.5 and m,h a = 70GeV. Here the coupling 
is restricted in the region between solid lines. As shown in the figure, Cntu * s a ^ ^he or der of 
the magnitude 10 2 and take smaller values if m^ takes larger values compared to rrih in the 
theory. If neutral Higgs masses are almost degenerate, namely m,A = m^, should be 

very large since the contributions of h and A have opposite signs and the same functional 
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dependences. Therefore to be able control the numerical value of the coupling £^ r/i , the masses 
should not be degenerate. 

Now, we use the upper limit of the BR of the process // — > e'y to restrict the Yukawa 
combination £Jv )TAt £jv )Te ( see ec i- (0))- Here we take into account only the internal r-lepton 
contribution. Using the restrictions for £^ T/J and CN,Tu,^N,Te we ^ n( ^ a constraint region for 
£^ re . In Fig. H we plot C^ Te as a function of for sinO^ = 0.5 and m/j = 70GeV. Here 
the coupling is restricted in the region between solid lines. As shown in the figure £jv Te * s a ^ 
the order of the magnitude 10 -4 and becomes smaller if the mass of the particle Aq is larger 
compared to m ho in the theory, similar to the coupling ^§^ T . Fig. [| represente the effect of 
possible mixing between CP even neutral Higgs bosons ho and H Q on the coupling £^ re . When 
the mixing parameter sina changes in the range — 0.1, the upper limit of £® re is affected 
at the order of the magnitude 3% for the fixed values of sinO^ = 0.5 , m ho = 70GeV and 
tuaq = 80 GeV . The lower limit of the coupling is almost non-sensitive to the mixing as shown 
in the figure. 

At this stage we are ready to predict the electron EDM, d e . Using eq. (Q) and the restriction 
for £^ re , d e is plotted with respect to sin^ in Fig. |5| for fixed values of sinO^ = 0.5, m ho = 
70 GeV and rriA = 80 GeV. For the intermediate values of sin6 e , d e is at the order of the 
magnitude 10~ 32 which lies in the experimental restriction region [|TJ| 

d e = 1.8 ±1.2 ±1.0 x 10~ 27 e- cm (18) 

We also present m ho dependence of d e for fixed values of sin9 e = 0.5 and rriA = 80Ge^ in 
Fig. ^. This figure shows that d e is bounded by the solid lines and decreases with increasing 
values of m^ . For completeness, we study h -H mixing effect on d e . Fig. [7| shows that d e is 
not sensitive to mixing of two CP-even bosons in the given range of sin a. 

Finally, we would like to predict the behaviour of d T with respect to BR(t — > fij). Since 
we take into account internal r and \i leptons in the calculation of d T , it has a functional 
dependence (see eq. (|J)) 

d T = h + h\H TT \ 2 , (19) 

where fi contains the coupling ^ T „ and its complex conjugate, which can be fixed by the 
constraint plotted in Fig. |2|. On the other hand BR(t — > ^7) has in the form 

BR(t - py) = 9i\Zn,t/ \Zn,tt\ 2 , (20) 
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where |£jv t J is again fixed by the constraint we have. As a result, we connect d T and 
BR(r — > /ry) by an expression in the form 

d T = fi + 92BR{r -> /ry) . (21) 

Here the lower limit of /1 is at the order of the magnitude of 10~ 15 for the fixed values of 
m-ho — 70GeV, rriA = SOGeV and sinO^ = 0.5. Therefore, even internal /i-lepton causes to 



exceed the upper limit of experimental value of d T ]T6 j 



d T = 3.1 x 10~ 16 e- cm , (22) 

In eq. (|2l|) , the part g2BR(r — > /ry) is due to the internal r-lepton contribution. Here #2 is 
at the order of the magnitude of 1CT 14 . Therefore it is not possible to get a strict bound for 



BR(t — ► /ry) using the eqs. (plf) and the experimental result eq. (|22|). 



In our work we choose the 2HDM type III and assume that only FCNC interactions exist at 

D 

N,Tfi' 



tree level. First we obtain constraints for the couplings £jv re and £jv T/ r The previous one is at 



the order of the magnitude of 10 4 and the last one is 10 2 . Numerical values of these couplings 
decrease when the mass difference of neutral Higgs bosons ho and Aq increases. Further, we 



predict the d e as 10 e — cm using the experimental result of d^ and the upper limit of 
BR(fi — > erf). This lies in the experimental restriction region of d e , however it is far from 
the present limits. Finally, we calculated the BR(t — > /ry) and d T and connect these physical 
quantities in the form of a linear expression (eq. (|2"ID). This expression does not provide a 
strict bound for BR{r — > /ry) and the constraint region for the coupling £^ Tr . 

In future, with the reliable experimental result of BR(r — > ^7) and obtaining more strict 
bounds for BR{ji — > e7), d e , and c? T it would be possible to test models and free parameters 
underconsideration more accurately. 
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Figure 1: One loop diagrams contribute to EDM of /-lepton and LFV interactions (if external 
leptons I have different flavors) due to the neutral Higgs bosons ho and A in the 2HDM. Dashed 
lines represent the electromagnetic field, h and A fields. 
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Figure 2: £j$ Tfl as a function of m^ for m ho = 70 GeV and sin 9^ = 0.5. Here the coupling is 
restricted in the region between solid lines. 
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Figure 3: The same as Fig. |] but for £ 
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Figure 4: £^ re as a function of sina for m ho = 70 GeV, m^ 
the coupling is restricted in the region between solid lines. 
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Figure 5: d e as a function of sin9 e for m ho 
in the region bounded by solid lines. 
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Figure 6: The same as Fig. |5] but d e as a function of m ho for sin9 e = 0.5. 
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Figure 7: d e as a function of sin a for mh = 70 GeV, rriA = 80 GeV and sin 6 e = 0.5. Here d e 
is restricted in the region bounded by solid lines. 
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